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Abstract 

We use a characterization of the fractional Laplacian as a Dirichlet 
to Neumann operator for an appropriate differential equation to study 
its obstacle problem in perforated domains. 

1 Introduction 

Given a smooth function ip : R" i-^^ R" and a subset of R", we consider 
v^{x) solution of the following obstacle problem: 

' v^{x) > ip{x) for X ^Tf, 

(-A)'*u^ > for X G R" (1) 
(-A)'u= = for X € R" \ and for x € if v^{x) > f{x). 

The operator (— A)^ denotes the fractional Laplace operator of order s, 
where s is a real number between and 1. It can be defined using Fourier 
transform, by ^((— A)'*/)(^) = \C\'^^f{0- In particular, ([1]) can be seen as 
the Euler-Lagrange equation for the minimization of the H norm 1 1 / 1 1^'' = 

1 1/(^)1^1* 1 1 with the constrain that / > (/? on Tg. We will see that this 
system of equations can also be stated as a boundary obstacle problem for 
elliptic degenerate equations. 

In ([T]), the domain R" is perforated and the obstacle ip{x) is viewed by 
v^{x) only on the subset Tg. A typical example of Tg is given by: 
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with a"^ <C e. The goal of this paper is to study the asymptotic behavior of 
as £ ^ 0. When is given by ([2|), the effective equation satisfied by the 
hmit of strongly depends on the radius a^: If is large enough, the limit 
turns out to be an obstacle problem with obstacle (p{x). On the other hand, 
if is small then the limiting problem is a simple elliptic equation without 
any obstacle condition. It is well known in the case of the regular Laplace 
operator (s = 1) that there is a critical size for a"^ for which interesting 
behavior arises. 

In the case of the regular Laplace operator, this problem was first studied 
for periodic Tg by L. Carbone and F. Colombini |CC80] and then in a more 
general framework by E. De Giorgi, G. Dal Maso and P. Longo [DGDMLSOj 
and G. Dal Maso and P. Longo [DML81j . G. Dal Maso |DM81j . Our main 
reference will be the papers of D. Cioranescu and F. Murat |CM82al ICM82b] . 
in which the case of a periodic distribution of balls is studied. More precisely, 
they prove that when s = 1 and if is given by ([2]) with = tq e"-^ ^ then 
the function v = lime_^o solves 

—Av — fi{v — if)- = 0, 

where /i is a real number (depending on ro) and W- = max(— 1(;,0). The 
obstacle condition thus disappears when e goes to zero, but it gives rise to 
a new term ^{v — ip)- in the equation. 

In [CM07] . we generalize this result (still with s = 1) to sets that are 
the union of small sets Se{k) C still periodically distributed, but with 
random sizes and shapes. More precisely, we introduce a probability space 
{Vl^J-,V) and we assume that for every w E and every e > we are given 
some subsets S'e(A;,u;) such that 

Se{k,Lo) C B,{ek). 

We then consider 

T,{u)= U Se{k,u:). 

The only assumptions necessary to generalize the result of D. Cioranescu 
and F. Murat |CM82aj - [CM82b] is that each set Ss{k, uj) is of capacity of 

n 

order e": cap{S£{k,uj)) = e^^{k,Ljj) (this is where the critical exponent e"-2 
comes from) and that the 7(/c, w) have some averaging properties (stationary 
ergodicity). 

In the present paper, we extend the result of |CM07] to the case of 
fractional Laplace operators s G (0, 1). We will show that under appropriate 
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assumptions on the size of the sets Ss{k, uj), the function v{x) = linig^o v'^{x) 
solves 

{-Ayv - ^i{v - (f)^ =0. 

In the particular case of sets of the form ([2]) , the critical size is now given 
by 

n 

(the critical exponent ^ "^s is related to the s-capacity of the sets ^^(A;,^;)). 



In the remainder of this section, we briefly motivate the problem and we 
introduce the extension problem for the fractional Laplace operators, which 
allows us to rewrite ([1]) as a boundary obstacle problem for a local (degen- 
erate) elliptic operator. The precise hypothesis on Ti.{uj) will be detailed in 
the following section in which the precise statement of the main theorem is 
also given. The remainder of the paper is devoted to the proof of our main 
statement. 



1.1 A semipermeable membrane problem. 

When s = 1/2, ([T]) naturally arises as a boundary obstacle problem for the 
regular Laplace operator (also know as Signorini problem): We consider the 
following problem set in the upper-half space M""^^ = {(x, y) G x M ; y > 
0}: 

-Au{x, y) = for {x, y) G W];^^ 
u(x, 0) > ip(x) for X £ 
dyu{x, 0) < for X G 
, dyu{x, 0) = for X G M" \ and for x G if u{x, 0) > ip{x) 

with the boundary condition 

lim u{x, y) = 0. 

y-^oo 

It is then well-known that v{x) = u{x, 0) is solution of ([1]) with s = 1/2 (see 
[SilOT] and [CSS07| for details). 

It can be of interest to state equation ([3]) in a bounded domain D C M"^^: 
Introducing 

S = D n {y = 0} and T = dD D {y > 0}, 
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we can consider the following boundary obstacle problem: 
—Au{x, y) = for (x, y) £ D 
u{x, 0) > ip{x) for X G S n T£ 



(4) 

dyu{x, 0) < for X G S 

dyu{x, 0) = for X € S \ and for x G if n(x, 0) > ip{x) 

with the boundary condition 

u{x, y) = g{x, y) for (x, y) G T. 

Equation ^ arises, for instance, in the modeling of diffusion through 
semi-permeable membranes (such as the membrane of a cell) : The membrane 
is modeled by the surface S. The outside concentration of molecules is given 
by ip{x), and the transport of molecules through the membrane and in the 
direction of the concentration gradient is possible only across some given 
channels (represented by the set T^) and only from the outside of the cell 
{{y < 0}) toward the inside of the cell D. At equilibrium, the concentration 
inside the cell is then given by the solution u{x,y) of 

1.2 An extension problem for fractional obstacle problems 

Following L. Caffarelli, S. Salsa and L. Silvestre jCSSOTj . we can actually 
rewrite ([1]) as a boundary obstacle problem for all fractional powers s G 
(0,1). We rely for this on the following extension formula established by 
L. Caffarelli and L. Silvestre [CS06] : For a given function /(x) defined in 
R", if we define u{x, y) by 

f -div (y"Vti) = for (x, y) G M"+^ 

\ (5) 
1 u{x,0) = /(x) for X G M", 



then 



(-A)^/(x) = Ihn y^dyu{x,y) 



with 

s = (1 - a)/2. 

We can thus rewrite the fractional obstacle problem ([1]) as follows: 
-div (y'^'Vu^) = for (x, y) G M++^ 
u^(x,0) > (p{x) for X G 

lim y"'dyu'^{x, y) < for x G M" 

lim y^dyu^x, y) = for x G M" \ and x G n {u^ > if} 
, y^o 



(6) 
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where a = 1 — 2s (note that a € ( — 1,1)). Our main result will concern 
problems such as ([6]) with possibly bounded domain D instead of M"^^. 

In the sequel, the theory of degenerate elliptic equations in weighted 
Sobolev spaces will play an important role. We refer to [FKS82| for many 
results that will be used in this paper. 

1.3 Variational formulation 



The system of equations ([6]) can also be written as a minimization problem. 
For a given opei 
space with W' 
space. We have 



For a given open subset D of R" , we denote by L'^{D, \y\'') the weighted 
space with weight \y\'^ and by VF^'^(D, |y|") the corresponding Sobolev's 



Fllvi/i.2(D,|j/|") 



/ \y\°-\u\^ dx dy + I \y\°-\Vu\^ dx dy. 
Jd Jd 



We then introduce the energy functional: 

^(u) = [ l\ynVu\'^dxdy 
Jd ^ 



and the set 

K^ = {v^ H^^'2(L>, lyl") ; f (x, 0) > ip{x) for x € T^X^) , v = g onV}. 

It is readily seen that © is the Euler-Lagrange equation associated to the 
minimization problem: 

^{n^)= ini^iv), u'GK,. (7) 

(Note that since is closed, convex and not empty, ([7]) has a unique 
solution S K^). 

Finally, we notice (see |CS06j ) that if u{x, y) is the extension of a function 
f{x) as in ([5]), then 

|a|V7..|2 j„, _ / \d\2s 



\ynVu\Uxdy= / \inf{OVdi = \\f\\^s 

In particular, the minimization problem ([7]) is equivalent to the variational 
formulation of problem ([1]). 
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In this paper, we study the asymptotic behavior of the solutions of ^ 
for any open subset D of M^"^^. The assumptions and the main result are 
made precise in the next section. The proof of the main theorem, which is 
details in Section [3l relies on the construction of an appropriate corrector. 
This construction is detailed in Sections H] and O 

2 Assumptions and Main result 

2.1 The set T, 

We consider a probability space {Q,J^,V). For all w € 17, the set T^ito) is 
given by: 

where the sets Se^kjUj) C M" satisfy the following assumptions: 

Assumption 1: For all k G anduj € 17, there exists 'j{k,uj) (independent 
of e) such that 

caps{Se{k,uj)) = e"7(A;,u;), 

where capg{A) denotes the s-capacity of subset A o/M""^^ (defined below). 
Moreover, we assume that 

Se{k, uj) C B^^^_^ (ek) for all /c G a.e. u e n, (8) 

for some large constant M , and that there exists a constant 7 > such that 

lik,uj) < 7 for all k G and a.e. uj £ 0,. (9) 

This first assumption defines the critical size of the set . It will guar- 
antee that capg(T(r) remains finite as e goes to zero. A natural definition for 
s-capacity of a subset A of M" is the following: 

cap,(^) = inf I / / G H^(Ml. /(^) > 1 for x e a] . 

Using the extension problem for the fractional Laplce operator (see |CS06j 
for details), an equivalent definition (up to a multiplicative constant) is given 
by 

cap,(^)=inf| / y'^lV/ipdxdy; /iGVFo'^(M++^lyr), /i(x,0) > 1, X G a1. 

UR"+i J 
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We will use this second definition in this paper. If is a n-dimensional 
ball, then its s-capacity in M"+^ is given by 

CaPs[J:!j. ) — Cn+l-ar — Cn+2sr 

for some constant c^. Assumption 1 is thus satisfied in particular if the sets 
S^{k,Ljj) are balls centered on eZ" with radius r{k,(jj)£n=^ . 

Assumption 2: The process 7 : Z'^ x i— >■ [0, 00) is stationary ergodic: 
There exists a family of measure-preserving transformations : ^ ^ Q, 
satisfying 

7(fe + k' , oj) = 7(fc, Tfc/w) for all k, k' € Z" and w G O, 

and such that if A C Q and t^A = A for all k G Z", then P{A) = or 
P{A) = 1 (the only invariant set of positive measure is the whole set). 

This second assumption is necessary to ensure that some averaging pro- 
cess occur as e goes to zero (the hypothesis of stationarity is the most general 
extension of the notions of periodicity and almost periodicity for a function 
to have some self-averaging behavior). 



2.2 Main result 

We are now ready to state our main result: 

Theorem 2.1 Let D be a open subset o/M""^^ (n>2), denote 

S = D n {y = 0}, V = ^DC^{y>Q} 

and let T^{(jj) he a subset o/S satisfying Assumptions 1 and 2 above. 

There exists a constant ao > such that for any (pix,y) e C^'^{D) the 
solution {x , y , Lo) of 

minj^y y^lVvf dxdy; v € Wq'^{D, |?/|"), v{x,0) > f{x,0) for x & T^{oj)^ 

converges t^^'^(Z?, \y\°')-weak and almost surely lo E to a function u{x, y) 
solution of the following minimization problem 

mm!^^Jjj''\Vv\^dxdy + ^J^aoiv - ^)lix,0)dx; v G W^'\d, {y^^ , 
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where w_ = max(0, — w). 

If, moreover, there exists 7 > such that j{k, cj) > 7 for all k ^T/^ and 
a.e. u; € fi, then oq > 0. 

In particular the function n(x, y) solves 

{-div (y^Vu) = for (x, y) ^ D 

lim y°'dyu{x, y) = ao(u — 0) for x G S 

u{x, y) = for (x, y) € F 

Remark 2.2 When D is a bounded subset of W^'^ , the condition u G 
Wg^'^(L>, lyl'*) cou/d easi/y be replaced by 

u G W''''^{D, \y\"-) , u{x, y) = g{x, y) for x G 9-D n {y > 0} 

for some function g{x,y) G L°°{dD n {y > 0}). 

We stated Theorem 12. II in its most general form. It contains the semiper- 
meable membrane problem, as well as our original problem ([1]) with the 
fractional operator. More precisely, if we have D = ]R!f.+i and if we consider 
the trace v{x) = ti(x,0) in Theorem 12.11 we get: 

Corollary 2.3 Let be a subset ofW^ (n > 2) satisfying Assumptions 1 
and 2 above. There exists ao ^ such that for any ip{x) G C^'^{W^), the 
solution f^(x,a;) of (OP converges, as e ^ 0, H^{D)-weak and almost surely 
to a function U(x) solution of 

{-AYv - aoiv - if)- = 0. 

As in Cioranescu - Murat |CM82a[ [CM82b] and Cafarelli-Mellet |CM07j . 
the proof of Theorem 12.11 relies on the construction of an appropriate cor- 
rector. More precisely, we use the following result: 

Proposition 2.4 Let T^[lo) be a subset of M" satisfying Assumptions 1 
and 2 above. There exists a non-negative constant oq such that for ev- 
ery bounded subset D ofMJ^^, there is a function WQ{x,y,uj) defined in D 
and satisfying 

w'ix, 0) = 1 /or X G Te{uj) f] {D n {y = 0}) (10) 
lk1lL-(D)<C (11) 
— >0 W^''^{D, \y\'')-weak a.s. u e n (12) 
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and 



For all sequences v^{x,y,uj) satisfying: 

' v''{x,0) > forx£ Te{uj) n S 



11^ 



^"loo(d) < C 



v'^ — > V in W^''^{D, lyl"") — weak, a.s. 
and for any (f) € T^{D) such that 4> > 0, we have: 



(13) 



lim / y"-Vw'^ ■ Vv'^(j)dxdy > 



with equality if v^{x,0) = for x E n S. 

The proof of Proposition 12.41 will occupy most of this paper. We stress 
the fact that Assumptions 1 and 2 are sufficient but by no mean necessary 
to the proof of this Proposition. Any set T^iuj) such that Proposition 12.41 
holds would be admissible for Theorem 12.11 

The condition (I13[) may seem rather obscure and the next Lemma will 
suggest a nicer (but stronger) condition to replace it. However (jlSp is the 
condition that appears naturally in the proof of Theorem 12.11 

Lemma 2.5 Let D be a bounded subset ofMJ^^, and assume that sat- 
isfies 

— div (y^-Vw^) = 
w%x,0) = 1 

limy^o y''dyW%x,y) = ao 



liTay^oy"-dyW'^{x,y) < 



for {x, y) e D 
forx e Te{uj)f\Y. 
for X E n S 
for X G Tg n S 



(14) 



together with ( li jj) and ^2\) . Then ^3\) holds. 

This lemma also gives an indication of how to construct w'^{x,y,uj): We will 
look for a constant oq such that the solution of (jl4p converges to zero in 
W^'^{D, |y|'^)-weak. 

Proof: Let v" G L°°{D) D W^^'^{D, \y\'^) be such that v%x,0) > on 
Tjr n S and let (p he a smooth test function with compact support in D. 
Then, we have: 

= / diY{y''Vw'')(j)v'' dxdy 
Jd 

= - y"-<pVw'' ■ Vv" dxdy - / y"V(/> • \7w%'' dx dy 
Jd Jd 

- / lim{y"'dyw'^) v'^ (pdx — / lim{y"'dyw'^)v'^ (pdx. 
Jt.\t, Jt, y-*^ 
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Since lim.y^Qy^'dyW^ {x,y) < and v^{x,0) > for x € Tg, we deduce: 

y''(j)Vw'^ ■Vv'^ dxdy > - y"'V(j)Vw'^v'^ dx dy - / aov'^cpdx. 

D Jd Js^ 



> 



D 



y°-V<l)Vw^v^ dxdy- j a^v^ 4> dx.{lb) 



with equality if v^{x, 0) = for x E T^. In order to pass to the hmit in (|15p . 
we note that we have the following convergences: 

— >Q VF^'2(D, lyl'^j-weak a.s. to ^Q, 

and 

— > V L^{D, [y|")-strong a.s. w G 0. 



Hence the first term in the right hand side of (jl5p goes to zero. Moreover 
we have 

v^{-,^) — > u(-,0) i?''(S)-weak and L2(S)-strong a.s. uj eVt, 
so gives 

lim / y'^cj)S/w'^ ■ Vv"^ dxdy > — UQVcpdx. 
with equality if v^[x, 0) = for x G Tg. □ 
2.3 Related problems 

Before turning to the proof of Theorem 12.11 we briefly mention other results 
that follow from Proposition 12.41 If we consider energy functionals of the 
form 

J'{v) = - / y''\Vu\^ dxdy + / uhdx 
2 Jd Jt. 

for some h G then a proof similar to that of Theorem 12.11 shows that 

the homogenization of the following equation 

' v^{x) > ip{x) for X G Tg 
{-AYv^ > h{x) for X G M" 
L {-/^Yv" = h{x) for X G M" \ and on Te if f ^ > 99 

leads to 

{-/\Yv-aQ{v-^)- = h inM". 
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More interestingly, we can replace the constrain v'^ > ip on by a 
Dirichlet condition of the form = on Tg. This amounts to minimizing 
(v) in the convex set 



K, = {ve W^'^{Rl+\ \yn ; v{x,0) = for x G r,(aj)}, 
The corresponding Euler equation is 



£ 



( {-Ayv''{x)=h{x) forxeW^XT, 
1 v^{x) =0 for a; € ST^. 

We can then show that the solution v'^{x) converges to a function v{x) 
solution of 

{-AYv - ao{v -ip)=h in M". 



3 Proof of Theorem D 

In this section, we prove that Theorem 12.11 follows from Proposition I2.4[ 
For the sake of simplicity, we assume that D is a bounded domain in R"^^. 
This allows us to take the corrector w'^{x, y, to) given by Proposition 12.41 and 
corresponding to the domain D. When D is unbounded, we note that every 
integral involving w'^ is computed with a compactly supported test function 
(p. We can thus use, for each of them, the corrector corresponding to the 
domain supp (p. The final result is of course independent of . 

The maximum principle and the natural energy estimate easily give that 
is bounded in L°°{D) fi Wq''^{D, \y\°') almost surely. In particular, there 
exists a function u{x, u, u) such that 

1 2 

— > u Wq' {D, |y|") — weak a.e. lo e Q. 
In order to prove Theorem 12.11 we have to show that 

^a(u) = inf ^a{v) a.e. cj G $7 (16) 

where is the energy associated to the limiting problem, given by: 

/a{v) = ]^ j^y''\Vv\^dxdy + ]^ j^ao{u-ipf_dx. 
Equality p6|) will be a consequence of the following lemmas: 
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Lemma 3.1 For any test function (j) € V^D), we have 
lim / y°'\Vw^\^(j) dx dy = / aQ(j)dx. 

Lemma 3.2 Let be a bounded sequence in W^'^{D, \y\'^) n L°°{D). If 

in W^'^{D, \y\°-)-weak, 

then 

liminf ^(u^) > ^Q,(n). 



Proof of Theorem 12.11 

For any v € T>{D), we consider the function v + {v — ip)-W^ (note that this 
function satisfies the obstacle constrain). Its energy is given by: 



V + {v — (p)-w'^) 



D 



dx dy 



D 



|Vt;p + \V{v - ip)-\'^w'^ + \{v - ip)-\'^\Vuff 
{v — ip)^V{v — Lp)^w^Vw^ + VvSI{v — ip)-W^ 

+S/v{v — ip)-Vw^ dx dy. 

Lemma [3. II and the weak convergence of to in thus imphes 

\un^{v + {v - ip)-W^) = J^a{v). 

Morever, it is readily seen that the function v + {v — ip)-W^ belongs to 
K^. Since minimizes ^ on ifg, we deduce 



v + {v- > /{u"), 



and therefore 



Ja{v) > limsup ^{u^) for all v G V{D). 



On the other hand, Lemma 13.21 gives 



liminf ^(n^) > ^a(u) 
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and so 



for all V G V{D). 



Equality (jl6p follows by a density argument. □ 

Proof of Lemma I3.lt This first lemma is a straightforward conse- 
quence of (|13p: If we take = 1 — w^, we have v^{x,0) = for a; € T^, 
v^{x, y) bounded in L°^{D) and y) converges to 1 in W^''^{D, |y|")-weak, 
L'^{D, |y|'^)-strong, and almost surely u; € 0. We can thus use (jlSp . which 
implies 



for all (/> G V{D). □ 

Proof of Lemma l3.2t Following Cioranescu-Murat (see |CM82b] . Propo- 
sition 3.1), we evaluate the quantity 



for some test function z with compact support in D and then take the limit 
as e goes to zero. 

Using (fT2]) . we obtain: 




and so 






Lemma 13.11 yields 
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Property p^ . together with the facts that € L^{D) and (u^ — ip){x, 0) > 
for X G Te, imphes 

hm / y"-{z - (^)_ Vu^ • Vu;^ = hm / y"-{z - v3)_ V(u^ - (p) ■ Vw^ dx dy 

+ hm / y'^{z — Lp)^'\/(p -Vw^ dxdy 
> — / aQ{u — ip){z — ip)^ dx. 



It follows that for any test function z G T>{D) we have: 

liminf / | Vn^ ^ > 2 / y^Vu-Vzdxdy- / y"-\Vz\^dxdy 

Jd Jd 



—2 / ao{u — (p){z — (f)- dx 
— / ao(-z — v^)?, dx. 



We can now take a sequence z„ that converges to It strongly in W^''^{D, 
and such that z„(-,0) converges to u(-,0) strongly in Using the 

fact that (u — ip){u — (^)_ = — (u — (/?)?., we get 

liminf / | Vn*^ ^ dx dy > / y'^|Vup dx dy + / ao{u - ip)'^_ dx . 

which concludes the proof.D 



4 The auxiliary corrrector 

4.1 Notations and scheme of the proof 

We recall that 

= {{x, y) G X M ; y > 0}, 

and we fix a bounded domain D C M"^"*^. For any xq G M" and yo > 0, we 
introduce the following notation for the Euclidian balls: 

Br{xo,yo) = {(x,y) G ; (|x - xqP + |y - yol^)'/' < r} , 

B+(xo,0) = 5,(xo,0)n{y > 0}, 

B^{xo) = {x G M" ; |x - xoj < r} . 



14 



4.1.1 The fundamental solution 

We recall (see |CSS07] for details) that the function 



h{x,y) 



n-l + g 



with 



vrlr(Y) 



P + y 



solves 



— div (2/'^V/i)(x, y) = for i/ > 
\\my°-dyh{x,y) — > -6{x), 

where 5{x) denotes the Dirac distribution centered at in M". We also have 

dw{y''Vh) = -fin,aS{x,y) ^-™"+^ 



m 



where 5{x,y) denotes the Dirac distribution centered at in M""'"^ and for 
some constant ^n,a- 

4.1.2 An auxiliary corrector 

One of the key point in the proof of Proposition 12.41 is to see that away 
from e/c, the set Se{k,Lo) is equivalent to a (n+l)-dimensional ball. More 
precisely, we introduce the capacitary potential Lp\{x,y,uj) associated to the 
set Se{k^Lv). It is defined by the following minimization problem: 

inf I / y"|V(/Jpda;dy; ip € W^'^{Rl'^'^ ,\y\''), (/^(x, 0) > 1 Vx G Ss{k,uj)\ . 

It is readily seen that, almost surely in uj, (p\{x,y,uj) satisfies 

' -div {y^Vifl) = for (x, y) G Rl+^ 

V9|(x,0) = l ioi X e Se{k,uj) (17) 

limj^^o y"'dyipl{x, y) =0 for x ^ Seik, iv) 

and by definition of the capacity as seen in the introduction. Assumption 1 
yields 

y^\Vipl\''dxdy = e^j{k,uj). (18) 



/ 



n+l 



Moreover, we have the following lemma (the proof of which is presented in 
Appendix 
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Lemma 4.1 For any 6 > Q, there exists Rs such that 



2 



2 

< 5e'^^{k, oj) h{x — ek, y) 



for all (x, y) such that \{x — ek, y)\ > e^-'^+'^Rs and for all e > 0. 

Moreover, Rg depends only on the constant AI appearing in Assump- 
tion 1 (in particular, Rg is independent of k and uj). 

This Lemma will play a fundamental role in the proof of Proposition l2.4l (see 

n 

Section[5]). It suggests that at distance away from ek, the corrector 

should behave like the function 

2 

h%{x,y,^) ■■= e'^-i{k,uj) h{x -ek,y). 

For later use, we introduce the notation 



= £n-l + a _ 

The first step in the proof, and the main goal of this section is to construct 
a function w'^ that would be a good approximation of away from ek and 
that behaves like hf. at distance a^R from ek 

For that purpose, we introduce 

D, = D\\J 5+,^^)^, (ek) , and = S \ (ek) , 

where r{k,uj) is chosen in such a way that hf,{x,y) = 1 on dB^^^ ^y^(ek), 
i.e. 

X l/(n-l+a) 

r(A.,^)=P^7(A^,^) • (19) 

We will prove the following proposition: 

Proposition 4.2 There exist a non-negative real number oq (independent 
of the choice of D) and a function w^{x,y,u!) satisfying 

-div{y''Vw^) = for {x,y) G 

(20) 

lim y"'dyW^{x, y) = for x € He 

for almost all lo ^Vt, such that 

w'{x, y) = hl{x, y) + o(l) for {x, y) G B+^{ek) n De a.s. uj e n (21) 
Moreover, we have: 
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(^) II«^'IIl-(d,) ^ 

(a) \\w^\\i^20 ) — > as e ^ 0. 
(m)\\yw^\\L2 0^^ < c 

The goal of this section is to estabhsh Proposition 14.21 The main ad- 
vantage of over w'^ is that the former only depends on the capacity of 
Ss{k,uj). This explain why no assumptions are needed on the shape of 
Se{k,Lu). In the last section of the paper (Section [5]), we will see how to use 
both the functions (near ek) and the corrector w'^ (at distance a'^R of 
ek) in order to prove Proposition 12. 4[ 

4.1.3 Effective equation 

The main idea to prove Proposition 14.21 (and in particular (|21|) ) makes use 
of the fact that /i|,(x,y,a;) solves: 



— div (z/"V/i|)(x, y) = for {x,y) € 

lim y"'dyhl{x, y) = -e"-j{k, io)5{x - ek) for x £ M" 

with 

2 

7(A;,a;) = 'y{k,uj)- 



pn+l 



Proposition 14.21 will thus be a consequence of the following proposition: 
Proposition 4.3 There exists ao > such that the solution WQ{x,y,uj) of 
- div iy'^Vwff) = for (x, y) G Rl^'^ 

lim y^-dywl = ao - e'^^{k,uj)5{x - ek) for x eT, (22) 

fceZ"nE 

w§(x,0)=0 /orx€M"\S 
satisfies: 

w'o{x,y) = hUx,y)+oil) for (x,y) G B+^{ek)nD a.s. ugQ (23) 

This proposition is the main step in the proof of Proposition 14.21 and its 
proof will occupy most of section. 
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4.2 Proof of Proposition 14.31 

In order to prove Proposition 14.31 it is more convenient to work with the 
rescaled function 

voix, y, w) = e-^+^wliex, ey, w). (24) 
Equation (|22p then becomes: 

^ -div (y"Vv§) = for {x, y) G 

hm y"'dyVQ{x, y) = uq — > lik, uj)5{x — k) for x € e^^T^ 

fg(x,0) = for X G M" 

(25) 

and ([23]) is equivalent to 

VQ{x,y,u!) = hh{x,y,uj)+o{e~^'^'^) for {x,y) G -B^2(^^)'^^~"^-^ ^-S- <^ G 
where 

r(A;,w)"-i+'' 



hk{x,y) ■■= j{k,uj) h{x - k,y) 



|(x-fc)2+y2|(n-l+a)/2' 



1 _i_ 1 — g 

Note that hk = e~ on dB-^er(k,uj) with = e^-i+a . 

In order to find the critical ao for which the solution Vq has the appropri- 
ate behavior near the lattice points A; G Z", we follow the method developed 
by Caffarelli-Souganidis-Wang in |CSW05] and which was already the cor- 
ner stone in [CM07] : We introduce the following obstacle problem, for every 
open set A C M" and for every real number a G M: 

' v{x, 0) > for X G M" 

lim v{x, y) = for x G M" 

y-^oo 

-div {y^Vv') > for (x, y) G 

lim y^dyv{x, y) < a — > u>)S{x — k) for x ^ A. 

(26) 

We then define the smallest super-solution of the obstacle problem: 

Va,A{x-, y, uj) = inf {v{x, y) ; v solution of ([26]) }. (27) 
It is readily seen that the function Va,A satisfies 

-div (y^VvaA) = fo^^ i^, y) G M++^ 

lim y°'dyVa,A{x, y) = q— > l{k, uj)5{x — k) for x G A n {va,A > 0} 

(28) 
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and 

lim y°-dyVa,A{x, y) > for x G ^ n {va,A = 0}. (29) 



hk,a{x,y) = hk{x,y) - a I u-i —dx' (30) 



Remark 4.4 The function 

l) - 

iB!^(k) {\x - x'p + 2/2) — 
satisfies 

-diviy'^Vkk^a) = forxe M++^ 

lim y°'dyhk a{x, y) — > a — ^{k, uj)5{x — k) for x G B'^{k). 

It is radially symmetric around k and sup^^^^^i y^Q ha^k{x,y) < j-'n-i+a^ 
particular, the maximum principle and (2B\) implies that if Bf{k) C A, then: 

Va,A{x^ y, to) > ha,k{x, y, w) - r""^+" for {x, y) £ (k), a.s. (31) 

We now want to show that there exists a critical oq such that the fol- 
lowings hold: 

1. The solution of the obstacle problem i;a^A(a;, y, uj) behaves like ha,k{x, y, uj) 
near any point k £ A fl Z". 

2. The solution of (j25|) is not far from Va,A- 

For that purpose, we introduce the following quantity, which measures 
the size of the contact set along the boundary {y = 0}: 

Tna{A,u;) = \{x e A; Va,A{x,0,uj) = 0}| 

where \A\ denotes the Lebesgue measure of a set A in M"'. 

The starting point of the proof is the following lemma: 

Lemma 4.5 The random variable rfia is subadditive, and the process 

Tkm{A, u) = m{k + A, uj) 

has the same distribution for all A; E Z". 
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Proof of Lemma 14. 5t Assume that the finite family of sets (Ai)i^j is such 
that 

Aid A for all i G / 
Ai n = for all i / j 

\A-\Ji<.iAi\ = 

then Va,A is admissible for each j4j, and so I'a,Ai < ^ia,A- It follows that 
{U„,a(-,0,u;) =0}nAi C {U«,a,(-,0,u;) =0} 

and so 

m„(A,u;) = ^|{?;„,a(-,0,w) =0}n^i| 

which gives the subadditive property. Assumption 2 then yields 

Tkra{A,uS) = m(^,rfca;) 
which gives the last assertion of the lemma. □ 

Since ma(^,u;) < |^|, and thanks to the ergodicity of the transforma- 
tions Tfc, it follows from the subadditive ergodic theorem (see |DMM86] ) that 
for each a, there exists a constant such that 

ma(-Bt(0),u;) - 
hm ; — — = l\a) a.s., 

\Btm ^ ' 

where Bt{0) denotes the ball centered at the origin with radius t. Note that 
the limit exists and is the same if instead of Bt{0), we use cubes or balls 
centered at txo for some xq. 

If we scale back and consider the function 

w%{x,y,uj) = e^-" Va^B^_i{e-^xo){x/e,y/e,u}), in Bi{xq), 

we deduce _ 

^.^|{x;u;^(x,Oa;) = 0}i^ a.s. 

e^O \Bi\ ^ ' 

The next lemma summarizes the properties of l{a): 
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Lemma 4.6 

(i) i{a) is a nondecreasing functions of a. 

(a) If a < 0, thenl{a) = 0. Moreover, if the 'y{k,uj) are bounded from below, 
then i{a) = for a positive small enough (0 < a < C{'y) ). 

(Hi) If a is large enough (a > C{j)), then J(a) > 0. 

The proof of this Lemma is rather technical and of httle interest. It is 
presented in fuh details in Appendix |B1 Using Lemma 14.61 we can define 

ao = sup{a ; l{a) = 0}. 

We observe that ao is finite (Lemma 14.61 (iii)) and that ao is non negative 
(Lemma l4.6l (ii)). Moreover, ao is strictly positive if the 7(/c,u;) are bounded 
from below almost surely by a positive constant. 

We now fix a bounded subset A of M" and we denote by 

Vai^, y, = v^^^-iAix, y, t^) (32) 

the solutions of ([27l) corresponding to e~^A. We also introduce the rescaled 
function 

w%{x, y, io) = e^"" v%{x/e, y/e, w). 



In order to complete the proof of Proposition 14. 3^ we are first going to 
prove that w% satisfies inequality (p3]) , and then that the solution Wq of (p2]) 
behaves like IZJ^. 

We recall the definition of /i^ fc : 



hnk(X,V) = TZ Oi / TJ— ax , 

ct,n\ ^ a I n — l + g / „ n— l + g ' 

{\x-kY^y^) 2 J Bi{k) {\x - x'Y + y^) 2 

and we introduce the scaled function 

h'a,kix,y) := e"^""/iQ,fc(a;/e,y/e). 
Note that when (x,y) € dBl^^^^^ ^^{k), then 

7 / N -l+a / ^n+l+a , / 

K,k[.x,y)=e ^ -a dx 

Jb'-KQ) {\x - x'p + y2) 2 

_ l-a 

(we recall that a^ = 

We then have the following lemma: 
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Lemma 4.7 

(i) For every a and for every k (z 'Z"' f] A, we have 

vlix, y) > K^k{x, y) - r"-i+" for {x, y) G B+{k) a. s. 
(a) For every a > oq and every k £ X'^ A, we have 

vl,{x, y) < ha,k{x, y) + o(e"^+") for (x, y) G B^^^{k) a. s. 
We deduce: 
Corollary 4.8 

(i) For every a and every k "Z"' Ci A such that r{k,uj) > 0, we have 
vl{x,y) > e-^+^ + o(l) for (x,y) e dB+,^^^^.,{k) a.e.uj^n 
and so 

w%{x,y) > 1 + oie^-") for {x,y) G dB+^^^^^,{k) a.e. w G 
for all a. 

(a) For every a > oq and every k £ 'Z'"' f] A, we have 

vUx,y) < + for {x,y) G dB+j^^^^-,{k) a.e.coeQ 

and so 

wUx,y) < 1 + o(l) for {x,y) G (/c) a.e. LveQ 

Proof of Lemma 14. 7t 

(i) This is an immediate consequence of (|3ip . 

(ii) The proof of (ii) is more dehcate and is split in several steps. 
Preliminary: First of all since A is bounded, we have A C B^(xo) for some 
R. Without loss of generality, we can always assume that B'^{xo) = -Bf (0). 
If we consider 

vl{x,y,uj) =Va^^-iBn{x,y,uj), 
the solution of ([27|) corresponding to A = B"_i{0), it is readily seen that 

lfa{x,y,uj) < v^^{x,y,uj) for all (x, y) G M"+^ a.e. w G ri. 

It is thus enough to prove (ii) for v^. 

In the sequel, we will need the following consequence of Lemma 14.51 (see 
[CSWOSj for the proof): 
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Lemma 4.9 For any ball B^{xq) € i?"(0), the following limit holds, a.s. in 

|R(x,0,a;) =0}ni?;^^(£-ixi)| 
\B^\ = 



Step 1: We now start the proof: For any 6 > 0, we can cover S^ii by a 
finite number N (< Cd^^) of balls i?" = B^^_i{e^^Xi) with radius 6e~^ and 
center e~^Xj. Since a > ao, we have J{a) > 0. By Lemma IT9| we deduce 
that for every i, there exists such that if e < then 

|{t;^(x,0,u;) = 0} nSf I > a.s. w. 

In particular, if e < inf Sj, then v^{x'-,0) = for some x- in i?" a.s. l<j E fi. 

Introducing i?j = S^^-i (e~^Xi) the n + 1 dimensional ball with same 
radius and same center as Bf", we now have to show that remains small 
in each B^ as long as we stay away from the lattice points A; G Z". More 
precisely, we want to show that 

sup v'^{x,y) <CS^-''e~ 

Ukez"B+ik)\B+^{k) 



Step 2: Let r/(x) be a nonnegative function defined in M" such that < 
r]{x) < 1 for all x, r]{x) = 1 in -Bi/g and r/ = in R"\i?;^/4. We then consider 
the function 

U = vl-kx7] 

where -k^ indicates the convolution in M" with respect to the x-variable. The 
function u{x,y) is nonnegative on 2Bf and satisfies 

r div {y-Vu) = for (x, y) G 2B+ 

\ -C < limj,^o y^dyuix, y) < C for x G 2Bf ^ ' 

where C is a universal constant depending only on n, r and a. We deduce: 

Lemma 4.10 There exists a universal constant C such that 

supu < C7inf n + C5^-''e~^+'' . 
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Proof: We write u = U1 + U2 where ui and U2 are two functions solution of 
div {y"'S/ui) = in 2B^ and satisfying 

{lim y"' dyUi{x , y) = lim y"'dyu{x,y) for x G 2Sf, 
ui{x,y)=Q for (rc,y) €9(25+) n{y>0} 

and 

{lim y'^dyU2{x, y) = for x G 25", 
n2(x,y) = u{x,y) for (x,y) G 9(25+) n {y > 0}. 
The maximum principle and the fact that Bi has radius 5e'^ yield: 

\u,{x,y)\ < C{{26e-')'---y'--) 

for all {x,y) G 25^+. On the other hand, boundary Harnack inequality for 
degenerate elliptic equation (see |FKS82j ) implies 

supti2 < C inf U2- 

The Lemma follows easily.Q 

For the next step, we will need the following lemma: 
Lemma 4.11 If v satisfies 

diviy^Vv) = in B^{xo,0) 

and 

lim y^dyv{x,y) < a for x G B^^^xq), 

then 

2 f 1 

— - / \y\''v{x,y)dxdy <v{xo,0) + aC{n)r " 

cj„+ar'^+" Jb+(xo,o) 

where C{n) is a universal constant and tOn+a = Jbi{xq 0) \y\"' dxdy. 

Proof: The function w(x,y) = v(x,y) + a fnur ^ C'n+i+o^^ 

isfies 

div (y^-Vw) = and lim y^dyW < 0. 

y-+0 
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Proceeding as in |CS06| . we now reflect w about the plane {y = 0}. The 
function 

w{x, y) if y > 
w{x, —y) if y < 

is now defined in the whole space M""*""^ and it satisfies 



w{x,y) 



divdyl^VtU) < \nBr{xQ,Q). 
We can thus use the mean value formula (see |CSS07] ): 
1 



„,„ , .y\"'w{x,y)dxdy 

< IiJ(xo,0) 

< t(;(xo,0) 



< v{xq, 0) + a 

Je 



dx'. 



Since a > 0, we see that v < w and so 
2 



/ y°-v{x,y)dxdy < f \y\"-w{x,y) dx dy 

JB+(xn.O) ^n+ar J Br-(xn.O) 



^n+a'f'^^"' Jb+(xo,0) ^n+ar'"^" J Br {xo,0) 

Moreover, we have 

f Cn+l+a , I f Cn+l+a , n i-a 

/ 1 71 — n— dx = / — — r— dz = C [n + ajr , 

/„ , N ™/ ra— 1+a I , , ~ n— 1+a ^ ' ' 

hence the lemma. □ 



Step 3: We have Va{x[,0) = and \imy^Qy'^dyV^{x,y) < a for x G 
Bi/2{x'j). Lemma [4. Ill thus applies and yields: 

[ \y\%l{x,y)dxdy <C{v'^{xlO)+a) <C{a,n + a). (34) 

JB+^{x[,0) 

We want to deduce an upper bound on u in Bi. Since n > 0, we note that 

/■1/4 / \ /■1/4 

/ T°'u(x,T)dT > [ inf u] / T^dT. 

Jo Ve[o,i/4] J Jo 
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Then, using the definition of u (and the fact that ry(x) = outside By^{x)), 
we deduce: 

-1/4 



inf u < C inf / t"'u{x, t) dr 
t/M^o) ^ Jo 

/■1/4 r 

< Cinf / / TV^{^,T)d^dT 

^ Jo Jb",Ax) 



< C I T-V%{i,T)didT, 



Which, together with ([3^ yields: 

inf u < C(a, n). (35) 

Using Lemma I4.1UI we see that for every 5 and for e small enough, we have: 
supu < Cinfu + C5i-"e-i+" < C{a,n) + C5^-''e'^+'' < C5^-''e-^+'' . 

(36) 



Step 4: We now want to use (j36p to get an upper bound on f^. For that 
purpose, we note that lim^^o > in \ (~^k&"{k}-, and so a proof 

similar to that of Lemma 14.111 yields 

vlix^y) < Cn+a I |r|X(e,T)dedr (37) 

JB+^ix,y) 

for all (x,y) e Bi\ rikeZ"Bi/i{k). 

Inequality ([37|) and the definition of u{x,y) yield that for all (x,y) in 
Bi \ rikeZ" Bi/4{k), we have: 

ry+l/ 



<(x,y) < C„+, r / \r\%l{tr)d(dT 
ry+1/8 

< Cn+a / \t\"'u{x,t) dr 

Jv-l/S 



ly-1/8 

< C{n + a)\y\^^"' sup u. 



Inequality (j36p therefore implies 



sup <(x,y) < (38) 

(a;,s/)eUfcgznB+(fc)\B+4(fc) 
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Step 5: In order to complete the proof of the lemma, we only have to 
notice that since infg^^^^ y) > —Ca, ([38]) and the definition of 

imply 

y) < K,k{x, y) + C5^-"e-^+'^ in By^ik) 

for all G Z". □ 



This conclude the proof of Lemma 14.71 and we are now in position to 
complete the proof of Propositions 14.31 

Proof of Proposition 14.31 

For every a, we denote by the solution of the obstacle problem (j26p 
corresponding to ^4 = 

Va{x,y,^)= Va,e-^T.{x,y,^), 

and by w% the rescaled function: 

nf^{x,y,uj) = e^"" Va^^-iY;{x/e,y/e,uj). 

We recall that Wq is solution of 

-div iy^Vwf^) = for {x, y) G 

lim. y"'dyWQ{x,y) = ao — > j{k,uj)6{x — ek) for x G S 

wf^{x,0)=0 forxGM"\S 



In order to prove Proposition 14.31 we have to establish (j23p . This is done in 
two steps using the properties of the function w%: 

1. For every a > ao, we have div {y°'V{wQ — w%)) = for (x, y) G M!|:"'""^, 
lim y°'dy{wQ - w%) > ao - a on S and {wq — w%){x, 0) = on \ S. 

We deduce 

, , . _ , . f — a 

Wo{xo,yo) -w^{xo,yo) < / dx, 

Jt. |(xo - xY + y^\ 2 

and therefore 

l-g 1^ 

sup {wl{x,y) - w%{x,y)) < C|S|"+ip^+' \a - ao| 
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with 

PT. = inf{p; E C Bp}. 
In particular, we thus have 

wl{x,y) < + O(a-ao) for (x,y) G IR!^+\ 

and Lemma 14.71 (ii) (since a > oq) yields: 

Wq{x, y) < h% ^,{x, y) + 0{a - ao) + o{l) for (x, y) G B^i2{ek) a.s. 

(Note that this argument shows the continuity of with respect to 
a). 

2. Similarly, we observe that for a < oq, we have div {y"-S/{w^ ~'Wo)) = 
for (x, y) e M++\ « - wff){x, 0) = for x G R" \ S and 

- WQ){x,y) >a-ao- al{^^=o)nE for x G S. 

y—>u 

Proceeding as before, we deduce: 

sup < Cp^ \\J:\^{ao-a) 

+Ca\{wl{x,0) =0}nS|^ 

So Lemma 14.71 (i) yields 

w'oix, y) > h%^k{x, y) - o{e) - O{ao -a)- Ca\{wl{x, 0) = 0} n 
for all {x,y) G B^/2i^k). Finally, using the fact that 
lim|{<(x,0) =0}nS| = i{a)\T.\ = 

e— >0 

for all a < ao we easily deduce the first inequality in (f23l) . 

□ 



l-a 
1 
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4.3 Proof of Proposition 14.21 

In order to complete the proof of Proposition 14.21 we construct a corrector 
which is equal to 1 on the (n+l)-dimensional balls -S^^ ^^^^(eA;). More 

precisely, we recall that is a bounded subset of M""''^, and we introduce 



and 



r(fe,tij)a 



We then define a corrector w^{x, y, uj) which will satisfy all the conditions of 
Proposition 12.4^ with the set instead of T^. In particular, we will prove 
that behaves like M near the B^,, ^ ^iek). 
We consider the following obstacle problem: 

' div (y^V-u;) < for (x, y) G \ % 

lim y"'dyw(x, y) < Qn for x G Eg 

y-^o y V - u s ^ ^^^^ 

w{x,y) > 1 for {x,y) G 

^ w{x,0) =0 for X G M" \ S, 

and we define: 

w^{x, y, uj) = inf {w{x, y, w) ; w solution of ([39])} . 

It is readily seen that satisfies (pOj) . So in order to complete the proof 
of Proposition 14.21 we only have to show that w'^ is bounded uniformly in 
L°°{D) and that — > in W^^^{D, \y\"')-weak as e goes to zero. 

Strong convergence in L^{D, 

First of all, since = 1 = h'^^ ^(x, y) + o(l) on T^, (j23|) implies 

'WQ{x,y) — o{l) < w'^ {x,y,LLi) < Wq{x , y) + o{l) inD a.e. u £ Q, 
which in turn implies (using Proposition 14.31 again): 

hlf,{x,y)-o{l) < w%x,u;) < h^^ix , y) + 0(1) y{x,y) e B+^{ek). (40) 
In particular, we get: 

||'5J^||loo(M"+1) < C. 
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Moreover, a simple computation shows that 

JB,{ek)\B^e{ek) 

and it is readily seen that (j40p implies 

Hix,y)\<Ce'- + o{l) = oil) V(x,y)G (J dB,/^iek). 

We deduce: 

||5^"|li2(B,|j,|.) < Yl [ y^lhlkl^dxdy + oil) f \y\^dxdy 

and since T^jZ" n e^^S} < Ce"" for all n, we have: 

||tS^|li2(^,l,l„)<e + o(l) = o(l). (41) 

In particular 

w"^ — ^0 in L^(i:>, |y I") - strong. 

as e goes to zero. 
Bound in W^^'^{D, 

Using the definition if and an integration by parts, we get: 
/ ^ y''\Vw''\^ dxdy = / ^ y^V w"" ■ V {w'' - I) dx dy 

^ ^ i'^^ra y''wl{x,y)]{w^x,y) - I) da {x,y) 



= — ao / {w'^{x,0) — 1) dx 
The L°° bound thus yieds 

/ ^ y''\Vw'\'^dxdy < C7aolSel(l|ty"||L°° + 1) < C, 



which completes the proof. □ 
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5 Proof of proposition 12.41 

This section is devoted to the proof of the main proposition. We recall that 
the sets Sc{k,Lo) are subsets of M" with unspecified shapes and they satisfy 

cap^{Se{k,uj)) =e"7(A;,u;). 

Lemma |4. II gives the existence of a function such that 

' div (y"V(^) = for (x, y) G M!^+^ 

< ip{x,0) = 1 for X G Ss{k,uj) 

lim y"'dif{x,y) = for x ^ Se{k,uj) 

and we let oq and w'^{x,y,tLi) be given by Proposition I4.2[ 
We then have: 

1. For a given 5 > 0, Lemma [4.11 implies that for every G and G 57 
there exists a constant Rs{k,uj) such that 

M{x,u:) - hUx,y,u;)\ < 6hUx,y,u;) < 6^^^^ (42) 

in B^^efi^ \ B^sRgi^k) and for all e > 0. It is readily seen that for any 
R there exists ei (R) such that 

o=ii<e'"/4 foralle<ei. (43) 

for some a > 1. 

2. Inequality (j2ip in Proposition 14.21 implies that for given 6 and R, there 
exists £2(6, R) < £i{R) such that for all e < 62(6, R), we have 

\w'{x) - hUx,y,u;)\ < j^^^^ in B+^{ek). (44) 

Thanks to (03]), Inequality (01]) holds in particular in B^^^ j^\B^s j^(ek) . 

The corrector will be constructed by gluing together the functions 
(near the sets Ss{k)) and the function w'^ (away from the sets Ss{k)). The 
gluing has to be done very carefully so that the corrector satisfies all the 
properties listed in Proposition 12.41 For a given e, we define 5^ to be the 
smallest positive number such that (03]) and ([H]) hold with 6 = de and 
R = Rg^. From the remarks above, we see that 5s is well defined as soon 
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as £ is small enough (say smaller than £2{l,Ri)). Moreover, for any 5 > 0, 
there exists eo = £2{^-, R&) such that < S for all e < Sq.Ih particular 

lim 5r = 0. 

From now on, we write 

Re = Rs^ ■ 

In order to define w'^ , we introduce the cut-off function ri^{x,y) defined 
on D and such that 

= 1 for {x,y)eD\ [j B+,j^^{£k) 
r]e{x,y)=0 for(x,y)G \J K^R^^k). 

We can always choose r] in such a way that 

\Vris\ < C{a^Rsy^ and \At],\ < Cia'Re)-'^ 
for e B^a<ij^^{£k) \ B^sj^^{ek). We now set: 

w%x,y) = r]e{x,y)w^x,y) + {1 - Tje{x,y)) ^ ipl{x , y) 1 ^+ ^(^^^^{x , y) . 
It satisfies 

' iflix, y) for (x, y) G I) n 5+^^ (eA;) VA; G 
for(x,y)GZ)\ \J B+^^^^{ek). 

To simplify the notations in the sequel, we denote 

The properties of w'^ are summarized in the following lemma, which 
implies Proposition 12. 4t 

Lemma 5.1 The function satisfies the following properties: 
(i) w^{x,0) = 1 for X ^ Se and \ \w'^\\ioo(^d-) < C . 
(a) w'^ converges to zero in L^{D, \y\^)-strong as £ goes to zero. 
(Hi) w'^ is hounded in W^''^{D, \y\"'). 



w^{x,y,uj) = < 



32 



(iv) satisfies |I 
Proof: 

(i) Immediate consequence of the definition of since = 1 on Ss{k, uj) 
and and ipf. are bounded in L°°. 

(ii) Since S'e(A;,u;) C B1^^j^{£k), we have: 

for all such that \{x — ek,y)\ > a^M. Since 9?| < 1 in Ba^Mi^k), 

we get 



) 



+C E / brfe"^^^^ /i(x-eA;)j 

+C E e2n(^^e^^n+l-2(n-l+a)^n-l+a 

Using (j43p and the definition of a^, we deduce: 

(-) 2n — an 

Estimate (|4ip thus implies 

lk^llL2(D,|s,"|) < ll^'^llL2(D,[y'^|) + 11(1 - ??e)<^^llL2(D,|y''|) = 

and therefore 

— > L2(D, |y|")-strong. 
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(iii) Next, we want to show that is bounded in W^''^{D, \y\'^). First, 
we note that outside UfcgX"-B£/2(£^) we have Vw^ = S/w^ which is 
bounded in W'^''^{D, Next, we see that in 5^/2(^^)1 we have: 

Vw' = Vri,{w' - ifl) + 7],Vw' + (1 - r?,)V(/7| (45) 

Since and (f^ are both bounded in W^''^{D, |y|") (thanks to (fTHIl ). 
we see that in order to show that Vw^ is bounded in L'^{D, we 
only have to show that 



/ y'^lVveiw' -ip')\^dxdy<C. 
Jd 



For that purpose, we notice that ()42p and (j44p yield 

W -^k\< C—^^ in B2R,a-{£k) \ BR,a^(ek), 
and so, using the definition of r]e{x,y), we deduce: 
y^Y^rieiw" -ip'')\^dxdy 

< I y^\Vr^e{w'-^l)\^dx 

i ,-.r„nr-.^ J Bon^„£ (ek)\Bi}^„e (ek) 



D 



p2(n-l+a) 

where we used the fact that we can always assume that 5^ < 1 and 
i?e > 1. For latter use, we note that we actually proved 

/ y''\Vr]e{w^ - ip^)\^ dxdy — ^ when e ^ 0. (46) 

JD 



(iv) It remains to show that (jlSh holds. We only show the inequality (the 
equality follows easily). Let be a sequence of functions satisfying: 

' v^{x,Q) > for X G 

< \\v^\\l^{d) < C 

— >v in ^^^'^(1), |y|'') — weak. 
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Then for any € V^D), we have: 
y"-Vw^ ■ Vv^cjidxdy 

D 

= — y'^Vr]^ ■ Vv^{w'^ — ip'^)(j)dx dy — / y^'Vw'^ ■ Vv'^tp-q^ dx dy 
Jd Jd 

- / y"V(/?^ • \7v''(f) (1 - r]s) dx dy 
Jd 

= - y'^^Ve ■ "^v^iw^ - Lp'')<j)dxdy 
Jd 

+ / {Yvca y"" dyW^)v'^ (j) rji; dx + / {\mi y"" dyLp'^)v'^ (f) {1 — rjf,) dx 

Jt. Je 

+ / y'^Vw^ ■ V{(t)7]e)v^ dxdy+ / y^Vip^ ■ V(0(1 - dx dy 

Jd Jd 

where we used the fact that div (y'^'S/w^) = on supp ij^ and div (y'^'Vip^) = 
on supp (1 — 7?^). The first term goes to zero thanks to (j46p and the 
weak convergence of Vw^ in L^{D, and the boundary terms sat- 
isfy 

lim / {\\m.y°'dyW^)v^(l)r]^dx = lim / aQV^(l)r]^dx 



hrn / a(jV(f>r]dx 



and 



hm / {lim y"" dy(p^)v'^ — ri^)dx = lim / {lim y"'dy(f^)v^ — r]^) dx 

< 0. 

Finally, the last two terms can be rewritten as: 

y'^Ww^ ■V{(j)r]e)v^ dxdy + / y'^V ip^ ■ V {4>{l - rie))v'' dx dy 
D Jd 



y''V{w^ - if^) ■ (V%) v'^cp dx dy 

+ / y''v^r]eVw'' ■W4>dxdy + / y"- v%l - rje) W if'' ■ V (j) dx dy 
Jd Jd 
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Using the weak convergence of Vw^ and Vip^ to zero, we see that in 
order to prove (|13p . it only remains to prove that 

/ y"-V{w'^ - ip'^) ■ (Vr?e) v'^(j)dxdy — > when e ^ 0. 
Jd 

Since v'^ is bounded in L°°, it is enough to show that 

/ |?/|''|V(w^ - 99^)1 \Vrie\dxdy — > when e ^ 0. 



For that purpose, we recall that 

Re 



-ipl\< in \ B+^^ , 



and 

div (y-V(5J- - ^D) = for (x, y) G i?+ \ 5+^^/^ 

lim - y) = qq for x G SJ^.^ \ B^,^ ,^. 

J/— >u ^' 

In particular, interior gradient estimates (see |CSS07] ) implies 

in B^^^R, \ B^R,- We deduce: 

\ynViw' -ip')\\Vve\dxdy 
< I\ + \ynV{w'-ipn\\Vrj,\dxdy 

+ C{a'Re)-^~'' I ^ \y\''dxdy 



Hp 



36 



Using (|13|) and the definition of a^, we deduce: 

/ \ynV{iS' - ip')\ |V%| dxdy < C5e + Ce''^''-^l 
Jd 

which concludes the proof since cr > 1 and hm^-^o Se = 0. 

□ 

Acknoledgment: L. Caffarelh was partially supported by NSF grant DMS- 
0140338. A. Mellet was partially supported by NSERC discovery grant 
315596. 
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A Proof of Lemma 14.11 

We now turn to the proof of Lemma 14. 1[ We take k = and we recall that 
ipQ is the capacity potential associated to Se{0). It satisfies (fT7|) and ([18]). 
We then introduce the function 

G{x, ^, y, r) = h{x - -t) + h{x -^,y + T) 

which satisfies 

div^,^(|T|"V5,rG) = -IIn,aSix - ^, y - t) - Hn,aS{x -S.,y + T) 

and 

lim T''drG{x,^,y,T) = 

T-+0 

for all X, and y. If y > 0, we deduce that for any function (p{x, y), we have: 
/ T''W^,rG{x, t y, r)Vi,Mt r) dr. 

Jt>0 

= - [ div {r'^V^^rGix, e, y, r))^(e, r) d^ dr 

Jt>0 

- [ \imT''drG{x,^,y,TMC,0)dC 

= Hn,a^{x,y). 

Moreover, if ifQ{x,y) is the capacity potential associated to ^'^(O), then (fT7|) 
yields 

/ T'^V^^rGix, y, T)V5,,^§(e, r) d^ dr. 

Jt>0 

= - [ G{x, e, y, r)div (r'^V^^rm, r)) d^ dr 

Jt>0 

- [ Gix,C,y,0)lhnT''drip'o{tr)dC 
= _2/" h{x-(,y)\imT-drm,T)dt 

Combining those two equalities, we get: 

l^n,a^k{x,y) = -"^ / - ^,y)lim T''dr(Po{^,T)d^. 
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Next, we note that (fTSj) yields, after integration by parts and using (fTT]) : 
and therefore 



Voix,y) e"7(0)/i(x,y) 



— [ [Hx - y) - Hx, y)] lim T^dr^lit t) d^. 



In order to conclude, we recall that ^'^(O) C i?A/a'E(0) and so we have 
1^1 < Ma^ in the previous integral. If is such that > Ra^ with 

R > 8M, we deduce that for all ^ € -^^(O), we have: 

\Kx - C,y) - Hx,y)\ < sup \v^^yh{x - C,y)\\C\ 



< 



< 



< 



< 



sup 



n — a 
2 



^^Hx,y) 



-h{x,y) 



We can thus write 



'Pf){x,y) e"7(0)/i(x,y) 



CM 2 , , , 
< Kx,y) 



< 



R f^n,a 

CM 2 

R f^n,a 



5.(0) 

e"7(0)/i(x,y) 



limr»a.<^§(C,r: 



where the right hand side is bounded by 6- — e'^j{0)h{x,y) if R is large 
enough. □ 
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B Proof of Lemma 14.6 



(i) For a given set A, it is readily seen from the definition of Va,A that if 
a' < a, then Va\A is admissible for the obstacle problem with a: It follows 
that 

Va,A < Va',A for any a, a' such that a < a 

and so a 1-^ ma{A,io) is nondecreasing. The result follows from the defini- 
tion of £{a). 

(ii) If a is negative, then we have 

lim y"'dyVa,tB{x, y) < for x € M"'. 

Since Va,tB{x,y) > for {x,y) € M"^^, we deduce 

Vc,tB{x,0)>0 forxGM". 
It follows that iriaitB, to) = for all t > 0, and so i{a) = for all a < 0. 
If r{k,uj) is bounded below: 

r{k, w) > r > for all k G Z", a.e. a; G 

then, we define 

j.n—1+a I- I' I 
Loix, v) = — TT Oi dx' — Co 

+ 2 Jbji(o) (|x - x'l^ + y^) 2 



with 

where e denote any unit vector in M". In particular, we have 

ip{x, 0) = if |x[ = 1, 

and, if a is small enough 

(/3(x,0) > if |x| < 1, 

and 

ip{x,y) < if |x| = 1, y > 
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(we note that <^ is the sum of a term which is decreasing with respect to \x\ 
and one which is increasing). Since ip satisfies 

Um y'^dyip{x, y) = a — ^6{x) >a — 7(0, u))5{x), 
y^o — 

for all X G -B"(0), we deduce 

Va,tB {x,0) > ip{x,0) > in (0) . 

Since we can do this in any ball Bi{k), we must have ma{tB''\ijj) = for 
all t > 0, and so J{a) = for all a small enough. 



(iii) We consider the function 

j^—l+a 

(\xY + y^) 2 



a 



- dx' + C, 



'By(o) [\x - x'Y + y 

where the constant C will be chosen later. It satisfies 

limy°'dyil){x,y) =a-^d{x) < a -^{0,uj)5{x) ^x e B^{0), 



ip{x,y) 



C when \x^ + — > 00. 



and we note that '0(x,O) is radially symmetric. Moreover, when a is such 
that 



a 



''"+.'+°, , dx' > f"-^+» 



B^(0) |ei-x'|"-i+« 



then 



■0(x,O) < C when \x\ = 1. 

Since div (y^V-;/') = for y > and Ymvy^Qy"-dy'ip{x,y) = for x ^ -^1^(0), 
the strong maximum principle and Hopf Lemma yield that the minimum of 
il){x,y) is reached for y = and x € -B"(0), and with an appropriate choice 
of the constant C, we can always assume that this minimum is 0: 

inf ibix,y) = inf ih(x,0) = 



Finally, if a is such that 
1 



a 



L 



sr(o) 



|ei _ ™/|n-l+a 



|ei — X 



Hn—l+a 



dx' > r"-i+«(4'*+i - 1) 



then ip(x,0) reaches its minimum when = Ra with < 1/4. 
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We now consider the function ip{x,y) defined by: 

J ^{x - k,y) for {x,y) G B^j^{k) 

We clearly have 

lim y°'dyip{x, y) < ao — ^{k, oj)6{x — k) for x G B?,Ak) 
j/-»0 ' 

and 

lim y^'dyifiix, y) < ao for x G M" \ Ujt/5^4(A;'). 

In order to prove that tp is a supersolution for the obstacle problem, we only 
have to check that 

ipix -k,y)= inf -^{x - k', y) for (x, y) G dBt^{k) 

k' ' 

or equivalently 

ipix, y) = inf ipix - k', y) for (x, y) G 55+4(0). 

It is readily seen that this amounts to showing that 

n,;/ ,1 , / ^ ^ r/r' > 4"~^r"~^+" 

for all k E {0} and all G (?i?^4(0). This inequality is obviously 

satisfied if a is large enough provided we can prove that 



/I »n K— 1+0 „ n— 1+0 

BfCo) (|x - XX + y^) 2 — k - x'\^ + y^) 2 



dx' > 



for all G \ {0} and all (x, y) G ^^+,4(0). This is equivalent to 

1 



f I dx' > [ 

whick holds for all (x,y) G dBy^{0) since |A;| > 1. 
By definition of Va,tB, we deduce that 

Va,tB{x) < (p{x,y) in tS" a.s. 
42 



dx' > 



In particular, this implies that Va,tB'^ vanishes in tB^ \ UkeZ'^Bi/2{k), and 
so 

m^itB^,iu) ( |Cl|-|^r/2l ^ _ UJn 

We conclude 

> 1 _ ^ > 0. 

□ 
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